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Abstract 

We have searched for r lepton decays, t~ Avr" and r~ An~ , that simultaneously violate 
the conservation of the lepton (L) and baryon (B) number using a data sample of 154 fb~^ collected 
with the Belle detector at the KEKB e~^e~ asymmetric-energy collider. No evidence for a signal 
was found in either of the decay modes and we set the following upper limits for the branching 
fractions: B{t~ Avr^) < 1.4 x lO"'^ and B{t~ Avr") < 0.72 x 10"^ at the 90% confidence 
level. This is the first search ever performed for these modes. 

PACS numbers: ll.30.Fs; 13.35.Dx; 14.60.Fg 
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INTRODUCTION 



While the Standard Model (SM) assumes both the baryon (B) and lepton number (L) 
conservation, in some extensions of the SM, such as supersymmetric 0] and superstring [2| 
models, B and L violation is expected while their difference B — L should be conserved. Ex- 
isting searches for baryon number violation are restricted to experiments with nucleons 
Limits from searches for proton decay imply that the r lepton decays with baryon number 
violation can be expected to have a probability well below the current experimental sensi- 
tivity 0]. However, high luminosity S-factories provide an opportunity to look for decays of 
the r lepton, D and B mesons that violate B and L with unprecedented sensitivity. Until 
now, searches for r lepton decays with a baryon final state involved a proton only j^. Upper 
limits for the branching fraction of these decays are in the range of 10~^ ~ 10"^. In a recent 
extension of the SM p, decays of r lepton as well as D and B mesons that violate B and 
L were considered using right-handed four-fermion couplings that conserve B — L, these 
suggested that r~ Air" decay involving the second and third generation transition might 
be interesting. 

We report here our search for Avr^ and r~ Avr^ decays with a data sample of 

154 fb~^ collected at the T{4S) resonance and 60 MeV below it with the Belle detector at 
the KEKB e~^e~ asymmetric-energy collider 0]. These r lepton decay modes have never 
been studied before [f] . 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere j^. 

Particle identification is very important in this measurement. It is based on the ratio of 
the energy deposited in the ECL to the momentum measured in the SVD and CDC, the 
shower shape in the ECL, the particle range in the KLM, the hit information from the ACC, 
dE/dX in the CDC and time-of-fiight from the TOF. To distinguish hadron species, we use 
likelihood ratios, for instance, V{p/7t) = Cp/{Cp + where is the likelihood for the 
detector response to the track with flavor hypothesis i. For lepton identification, we use 
likelihood ratios as electron probability V{e) Q and muon probability V{fi) ^3] determined 



by the responses of the appropriate subdetectors. 

For the Monte Carlo (MC) studies, the following pro gram s are used to generate back- 
ground events: KORALB/TAUOLA |ll| for r+r", QQ jlj for BB and continuum, BH- 
LUMI [i3 for Bhabha, KKMC Q for /i+zi" and AAFH [l^ for two-photon processes. 
Signal MC is generated by KORALB/TAUOLA. To take a range of possible r~ —>■ Ayr" 
and Avr^ decay mechanisms into account, we generate events under three different 

assumptions: uniform angular distribution in the r rest frame, V — A and V + A interactions. 
We initially assume all r decays to have a uniform angular distribution and take the other 
hypotheses into account later. The Belle detector response is simulated by a GEANT 3 
based program. Unless stated otherwise, all kinematic variables are calculated in the labo- 
ratory frame, while the ones calculated in the e~^e~ center-of-mass (CM) frame are indicated 



[t] Unless otherwise stated, charge conjugate decays are imphed throughout this paper. 
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by the superscript "CM" . 



DATA ANALYSIS 

We search for t^t~ events in which one r decays into Avr" and Avr" (signal side) and the 
other r (tag side) decays into one charged particle with any number of additional photons 
and neutrinos. Thus, for the B — L conserving modes the experimental signature is 

{t~ {p7i+) + TT"} + |r+ (a track)+ + {n^^'^ > 0) + X(missing)|, 

and for the B — L violating modes, it is 

{t- (pn-) + TT"} + |r+ (a track)+ + {n^^'^>0) + X(missing)|. 

Here charged tracks and photons are required to be reconstructed within the fiducial volume 
defined by —0.866 < cos 6* < 0.956, where 6 is the polar angle with respect to the direction 
opposite to the e"*" beam. Charged tracks should have momentum transverse to the beam 
Pt > 0.1 GeV/c and photons should have energies > 0.1 GeV. We denote the pion from 
r Att as tti and that from A — pvr as . We can distinguish between the B — L conserving 
and violating modes by the charge of these pions: the B — L conserving decay modes have 
an opposite sign combination on the tti and 712 charges, while the B — L violating modes 
have a same sign combination. 

We first demand that the four tracks have a zero net charge. The magnitude of the 
thrust ^3] is required to be larger than 0.9 to suppress the qq continuum background. The 
event should have a 1-3 prong configuration relative to the plane perpendicular to the thrust 
axis. We select A candidates via the pn~ decay channel based on the angular difference 
between the A flight direction and the direction pointing from the interaction point to the 
decay vertex (see Ref. 3| for more details). The proton from the A decay is identified by 



a condition V{p/7r) > 0.6, which has a 70% efficiency and 11% and 4% fake rate from K 
and TT, respectively. In order to avoid fake A candidates in which an e~^e~ pair from a 
photon conversion results in two tracks in the signal side, an electron veto {V{e) < 0.1) is 
imposed on the three tracks on the signal side. The reconstructed A candidate mass should 
be within ±5 MeV/c^ of the nominal A mass and its momentum, p^, is required to be 
larger than 1.75 GeV/c to reduce contributions from the generic t~^t~ and qq continuum 
background as shown in Fig. ^ A discrepancy between the data and MC observed in the A 
momentum distribution in the CM system shown in Fig. ^ (right) is apparently caused by 
the imperfect simulation of baryon spectra in the QQ generator. Since the final estimate of 
the background uses information from the sideband data, this discrepancy does not directly 
affect our results. 

The total visible energy in the CM frame E^^, defined as the sum of the energies of 
the A candidate, two other tracks and photon candidates, is required to satisfy 5.29 GeV 
< E^^ < 10.5 GeV (see Fig. 121(a)). To ensure that the missing particles are only neutrino(s) 
rather than photons or charged particles that fall outside the detector acceptance, we require 
the residual momentum vector Pmiss to have a magnitude above 0.4 GeV/c and to point into 
the fiducial volume of the detector: —0.866 < cos^miss < 0.956 (see Fig. |21 (b) and (c)). We 
select events with missing particles belonging to the tag side by requiring the opening angle 
between the residual momentum vector and tag-side track to be smaller than 90 degrees 
in the CM system (see Fig. El (d)). To suppress the continuum background, the following 
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FIG. 1: Reconstructed A candidate mass (left) and momentum (right) in the CM frame. The 
signal MC distributions are indicated by the filled histogram, all background including t+t~ and qq 
by the open histogram, and closed circles are data. While the signal MC distribution is normalized 
arbitrarily, the data and MC are normalized to the same luminosity. The selected area is indicated 
by the vertical lines with arrows. 



requirements are imposed: the number of the photon candidates on the signal and tag side, 
j^siG ^ ]^ g^^^ n^^^ < 2, respectively. Both the proton veto (V^p/n) < 0.6) and kaon veto 
{V{K/n) < 0.6) are applied to tti and the tag-side track. 

Finally, a condition is imposed on the relation between Pmiss and the mass squared of the 
missing particle, m^j^g. The latter is defined as -E^iss ~ Pmiss^ where -Emiss = -E'totai — -^^vis, 
-Etotai is the sum of the beam energies (11.5 GeV) and E^is is the sum of all visible energy. To 
further suppress background from generic r+r^ and continuum background, we require the 
following relation between the missing momentum Pmiss and m^j^g : Pmiss > 1.5[GeV~^c^] x 
m^jgg — 1.0[GeV/c] (see Fig. |2)). While this condition retains 89% of the signal events, 
81% of the generic t~^t~ and 77% of uds continuum background are removed. There is no 
significant difference between signal distributions for the B — L conserving and violation 
modes according to the MC. After applying the above requirements, 82 and 62 events are 
retained in the data for the B — L conserving and violating mode, respectively. 



RESULTS 

Signal candidates are examined in the two-dimensional space of the Att^/ Att" invariant 
mass, Minv, and the difference of their energy from the beam energy in the CM system, AE. 
A signal event should have Minv close to the r-lepton mass and A£' close to 0. The Minv 
and AE resolutions are parameterized from the MC distributions around the peak using 
an asymmetric Gaussian shape to account for initial state radiation with widths o^f^^l^^"" = 

4.6/4.0 MeV/c^ and aif'^"^ = 22/29 MeV, where the "high/low" superscript indicates the 
higher/lower side of the peak. 

We blind a region of iSciMin^ around the r mass in Minv and a region of —0.5 < AE < 0.5 
GeV so as not to bias our choice of selection criteria. Figure |3] shows scatter-plots for data 
and MC samples distributed over ±15a in the Minv — AE plane. For the B — L conserving 
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FIG. 2: Kinematical distributions used in the event selection after A mass and momentum selection: 
(a) the total visible energy in the CM frame; (b) the momentum of the missing particle; (c) the 
polar angle of the missing particle; (d) the opening angle between the missing particle and tag-side 
track in the CM frame. The signal MC distributions are indicated by the filled histograms, the 
total background including r+r~ and qq is shown by the open histogram, and closed circles are 
data. While the signal MC distribution is normalized arbitrarily, the data and MC are normalized 
to the same luminosity. The selected area is indicated by the vertical lines with arrows. 



mode, the number of data events and the total number of events expected from MC outside 
the blinded region (bounded by the vertical dotted line in Fig. |3] (a) and (b)) is 5 and 
8.5 ± 3.1 events, respectively, while that for the B — L violating mode is 6 and 5.1 ± 2.3 
events, respectively, indicating good agreement between the data and MC describing the 
background. The surviving background events are due to generic t^t~ decays (about 1/2) 
and uds continuum (about 1/2). The former events are dominated by the r — > 01(1260)2^,- 
decays, in which two of the three charged pions from the ai(1260) decay form a fake A 
candidate. The continuum background events have one true A that forms a signal candidate 
together with another track. 

As a signal region, we take an elliptically shaped area that contains 90% of MC signal 
events remaining after all the requirements in M\^^ — AE plane, as shown in Fig. |3] (a) and 
(b). It corresponds to a signal detection efficiency of 11.8 ± 0.1% for the B — L conserving 
and 11.7 ± 0.1% for the B — L violating modes. 
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FIG. 3: Pmiss vs rfi^jgg plots for (a) signal MC events, (b) data, (c) qq {uds) continuum MC events 
and T~^r~ MC events. The area indicated by the arrow is the selected region. 



We assume the background distribution to be fiat along the Mjnv axis, and then obtain 
the expected background in the ellipse to be 1.7±0.8 events for each of the two modes, using 
sideband regions inside the horizontal dotted lines and outside the vertical dotted lines, as 
shown in Fig.|3](a) and (b). We open the blinded region and find only one data event in the 
ellipse for the B — L conserving mode and no data events for the B — L violating mode (see 
Fig. m (a) and (b), respectively). Since no statistically significant excess of data over the 
expected background in the signal region is observed, we apply the frequentist approach to 
calculate an upper limit for the signal yield . The resulting limits for the signal yields at 
90% confidence level, sgo, are 2.76 and 1.30, respectively. The upper limits on the branching 
fraction before the inclusion of systematic uncertainties are then calculated as 

where A',-^ = 137 x 10^ and B{K J'vr") = 0.639 0]. The resulting values are B{t~ —>■ 
An') < 1.3 X 10"^ and B{t~ An') < 0.64 x 10"^. 
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FIG. 4: Scatter-plot of data in the Mjnv ~ Aii^ plane: (a) and (b) correspond to the ±15a area 
for the B — L conserving and violating modes, respectively. The 90% elliptical region shown by a 
solid curve in (a) and (b) is used for evaluating the signal yield. In (a) and (b), the vertical dashed 
lines denote the boundaries of the blind regions, while the regions inside the horizontal dotted lines 
and outside the vertical dashed lines are sidebands used to estimate the expected background in 
the elliptical region. Closed circles correspond to the data. The filled boxes show the MC signal 
distribution with arbitrary normalization. 



The sources of systematic uncertainty in the detection sensitivity 2eNT-^B{A —>■ p7r~) 
are listed in Table lU The uncertainty from A reconstruction is estimated to be 6.0% from 
the difference of the proper time distributions for data and MC simulation. The proton 
identification in A decay contributes another 3.0%, and the branching fraction B{A Ptt~) 
has an uncertainty of 0.8% The total uncertainty from the tracking reconstruction 
efficiency is 4.2%. Other sources of the systematic uncertainties are: the trigger efficiency 
(0.5%), selection criteria (4.0%), MC statistics (0.7%) and luminosity (1.4%). The trigger 
efficiency is estimated using the trigger simulation. Assuming no correlation between them, 
all these uncertainties are combined in quadrature to give a total of 9.1% in the detection 
sensitivity. An additional systematic uncertainty due to the estimation of the expected 
background includes those due to its statistical error and the background shape. By varying 
the assumptions about the background shape for the MC events, we checked that this effect 
is negligible compared to the background statistical error. 

While the angular distribution of r~ — * Avr" and r~ Avr" decay is initially assumed 
to be uniform in this analysis, it is sensitive to the lepton flavor violating interaction struc- 
ture [i^], similar to other analyses jlJ]. The spin correlation between the r lepton in the 
signal and that in tag side must be considered. A possible nonuniformity was taken into 
account by comparing the uniform case with those assuming V — A and V + A interactions, 
which result in the maximum possible variations. No statistically significant difference in 
the Minv - Ai? distribution or the efficiencies is found compared the case of the uniform 
distribution. Therefore, systematic uncertainties due to these effects are neglected in the 
upper limit evaluation. 
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A reconstruction 


6.0 


Froton identincation 


3.0 


Branching traction o[A — > pvr j 


0.8 


Tracking reconstruction efficiency 


4.2 


belection criteria 


4.0 


Trigger efficiency 


0.5 


Luminosity 


1.4 


MC statistics 


0.7 


Total 


9.1 



TABLE I: Sources of systematic uncertainties for the detection sensitivity 2£Nt-t-B{K pn ) (in 
%) 



The upper limits on the branching fractions at the 90% CL. taking into account sys- 
tematic uncertainties on both expected background and the detection sensitivity are then 
calculated using the POLE program 2^1, which integrates over Gaussian systematic errors, 
to be: 



B{t- 



An' 
An' 



< 1.4 X 10"' 

< 0.72 X 10" 



CONCLUSION 

In conclusion, we have searched for the decay modes that violate both lepton and baryon 
number conservation: r" Avr" [B — L conserving) and r~ Avr" {B — L violating) using 
data collected by the Belle detector at the KEKB e^e~ asymmetric-energy collider. We 
found no signal in the either mode. The following upper limits on the branching fractions 
were obtained: B{t- Att") < 1.4 x 10"^ and B{t- Avr") < 0.72 x 10"^ at the 90% 
confidence level. These are the first results ever reported for these modes. 
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